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Abstract 

Using squeezed vacuum states formalism of quantum optics, an approximate 
solution to the semiclassical Einstein equation is obtained in Bianchi type-I universe. 
The phenomena of nonclassical particle creation is also examined in the anisotropic 
\ background cosmology. 
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1. Introduction 

The present universe in its over all structures seems to be spatially homogeneous and 
isotropic but there are reasons to believe that it has not been so in all its evolution and 
that inhomogeneities and anisotropies might have played an important role in the early 
universe 1,2 . The isotropic model is adequate enough for the description of the later stages 
of evolution of the universe but this does not mean that the model equally suite for the de- 
scription of very early stages of the evolution of the universe, especially near the singularity 
I 2 . Also the most general solution of the problem of gravitational collapse turn to be locally 
anisotropic near the singularity 3_5 . Cosmological solutions of Einstein's general relativity 
are also known in which the expansion can be anisotropic at first, near the singularity, and 
later the expansion became isotropic. Also, to avoid postulating specific initial conditions 
as well as the existence of particle horizon in isotropic models, attempts have been made 
through the study of inhomogeneous and anisotropic models of the universe. Among the 
anisotropic cosmological models, Bianchi type -I universe is the simplest one. In this model 
the metric considered as spatially homogeneous and possibly anisotropic. In contrast to the 
Friedmann-Robertson- Walker (FRW) metric, Bianchi type-I metric has three scale factors 
which evolve differently in their respective direction. Therefor the expansion in this model 
could be considered as anisotropic expansion. Interests in such models have been received 
much attention 1 > 3 ~ 7 . Huang has considered the fate of symmetry in Bianchi type - I cos- 
mology using adiabatic approximation for massless field with arbitrary coupling to gravity 
8 . Futamase has studied the effective potential in Bianchi type - I cosmology 9 . Berkin has 
examined the effective potential in Bianchi type -I universe, for scalar field having arbitrary 
mass and coupling to gravity 10 . These studies show that Bianchi type -I cosmological model 
may be useful to understand the very early universe. Anisotropic models of the universe 
which become isotropic have been considered several times 11 . These motivate the study 
of anisotropic background cosmological models with scalar field possess the advantage of 
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FRW model, and to analysis the possibility of its approach to the isotropic model with 
accuracy required by observations 12 . With this aim Bianchi type -I universe has studied 
and the relevant equation for which extreme close in the form of isotropic case examined. 
Form anisotropic to isotropic transition a damping mechanism is required. One of the ef- 
ficient damping mechanisms could be due to the particle creation in anisotropic models. 
Therefore it would be useful to examine the particle creation in the anisotropic cosmolog- 
ical model with nonclassical inflaton, which could expect to produce sufficient particles to 
bring isotropy during the evolution process of the universe. 

In this paper we study a homogeneous massive scalar field (inflaton), minimally coupled 
to gravity, in a spatially homogeneous and anisotropic background metric. The inflaton 
under our consideration can be quantized and represent in squeezed vacuum states, hence 
an approximate solution to the semiclassical Einstein equation and the phenomenon of 
nonclassical particle creation can be examined in semiclassical theory of gravity. 

Throughout the paper, we follow the units c — G — % — 1. 

2. Inflaton in Bianchi type -I metric 

The Bianchi type- I model is an anisotropic generalization of the FRW model with Euclidean 
spatial geometry. The Bianchi type -I metric can be considered as spatially homogeneous 
and anisotropic and is given by: 

ds 2 = dt 2 - S 2 (t)d Xl 2 - S 2 2 (t)dx 2 2 - S 2 3 (t)dx 3 2 , (1) 

where Si(t), Sz(t) and 63 (t) are the scale factors in three spatial directions. Which are 
representing the size of the universe in their respective direction. The three scale factors 
Si(t), <S>2(£) and Ss(t) are determined via Einstein's equations. 

In the background metric (1), consider an inflaton, minimally coupled to gravity, satisfy the 
equation: 

(g^V,V u -m 2 ) ( f ) (x,t) = 0, (2) 

where V M is the covariant derivative. The Lagrangian density for the inflaton field, <fi, is 
given by: 

£ = -\v=g (g^d^d^ + m 2 2 ) . (3) 

Since the inflaton is homogeneous, ie.; <f)(x,t) = 4>(t), its classical equation of motion for the 
metric (1) can be written as 

^)+E(H|)^) + ^)=0. (4) 

In the present context (4) is the classical equation of motion for the inflaton for the metric 
(1)- 
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For the metric (1), the purely temporal component of the classical gravity is now the classical 
Einstein equation and can be written as: 

«Si(t)«S 2 (t) S 2 (t)S 3 (t) S^Ssjt) _ 8nT 00 

C /A C /A C /A C /A C /A C /A C /A ' v ' 



S l {t)S 2 {t) S 2 {t)S 3 {t) S 1 (t)S 3 (t) ^(^(t)^)' 



where 



— I I ; 
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Too = S l {t)S 2 {t)S 3 {t) ^ + m 2 ^ , (6) 



is the energy density of the inflaton. In cosmological context, the classical Einstein equation 
(5), means that the Hubble constants [h^ = are determined by the energy density of 

the dynamically evolving inflaton described by the classical equation of motion. 

In order to study the full quantum effects in a cosmological model, both metric and mat- 
ter fields are to be treated quantum mechanically. Since a consistent quantum theory of 
gravity is not available, in most of the cosmological models the background metric under 
consideration is taken as classical form and matter field treat as quantum mechanical. Such 
approximation of the Einstein equation is know as semiclassical approximation. 

In semiclassical theory Einstein equation takes the following form 

GV = 87r(T^>, (7) 

where G^ u = R^ u — ^g^R is the Einstein tensor and (T^) is the expectation value of the 
energy-momentum tensor for a matter field in a suitable quantum state under consideration. 
In (7) the quantum field can be represented by a scalar field, 0, and is governed by the time 
dependent Schrodinger equation 

i^>(<l>,t) = H m (<f>,t)*(<l>,t). (8) 

Using the canonical quantization procedure, the scalar field can be quantized by defining 
the momentum conjugate to (f>, as 

Thus the inflaton in the Bianchi type -I cosmological model, can be described by a time 
dependent harmonic oscillator with the Hamiltonian 

Hm = 2S 1 (t)S 2 (t)S 3 (t)^ + 2 Ht) - (10) 

The eigenstates of the Hamiltonian are the Fock states which can be constructed by anni- 
hilation and creation operators in the following manner. 

a(t) = (j)*(t)^ - S 1 {t)S 2 {t)S 3 (t)^*(t)l (11) 
a f (t) = 0(*)tt - Si(t)S2(t)S 3 (t)0(t)0. 
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Thus the Fock state of Hamiltonian is 



a) a | n, <p,t) = n \ n, (ft, t). (12) 
In the present context the semiclassical Einstein equation takes following form: 

s^t) s 2 (t) ^ s 2 (t) s 3 (t) ^ s.it) s 3 (t) s^s^s^ty h 1 J 

where H is given by (10). 

3. Particle Creation in the anisotropic universe 

Most of the cosmological models are based on a classical behaviour of the scalar field. 
Therefore, it is of interest to study the evolution of the system with the scalar field, which 
possess the nonclassical features. Recently squeezed states formalism of quantum optics 13 is 
found much useful to deal with many issues in cosmology 13-23 . Squeezed states are minimum 
uncertainty states and are obeying Heisenberg uncertainty principle. A squeezed state is 
generated by the action of squeezed operator on any coherent state and, in particular, on 
the vacuum state. Therefore a squeezed vacuum can be defined 13 

\Z) = Z(r,<p)\Q), (14) 

where the squeezing operator, 

Z(r, if) = exp (e~ itp a 2 - e^a f2 )) . (15) 

In (15) r is the squeezing parameter which determines the strength of squeezing and if is the 
squeezing angle which determines the distribution between the conjugate variables and they 
take values 0< r < oo and — 7r < tp < n. a and cr are respectively known as annihilation 
and creation operators. When the squeezing operator acts on annihilation and creation 
operators lead the following results 13 : 

Z^aZ = a cosh r — a'e l(p sinh r, (16) 
Z^a) Z = a) cosh r — ae~ llp sinh r. 

In the case of squeezed states, the variance of the quadrate components are not equal but 
one component of the noise is always squeezed with respect to another. 

Next, consider the Hamiltonian of the semiclassical Einstein equation (13), whose expecta- 
tion value can be computed in squeezed vacuum state by replacing the number state | n, (j), t) 
with | Z). Therefor using (14), (15) and (16) in (13), we obtain the semiclassical Einstein 
equation as: 



S 2 (t)S 3 (t) S^Ssjt) 
S l (t)S 2 (t) S 2 (t)S 3 (t) 5i(*)5 3 (t) 



(sinh 2 r + 1) (0*0 + m 2 0*0) + (17) 



Re jcoshr sinh re %Lp (J) 2 + m 2 2 )} 
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In the above equations, and (jf satisfy the boundary condition 

5!(t)5 2 (t)5 3 (t) (j>*(t)<f>(t) " = 



(18) 



To solve the semiclassical equations (17), transform the solution in the following form 

4>{t) = [Si(f)&(f)S s (f)]- 1/2 x(t) (19) 

Therefor (4) becomes 



m + 



rn 



&(f) 



The inflaton has a solution of the form 



x(t) = o. 



x(t) 



-^exp(-i/ 7 ( t )<tt), 



(20) 



(21) 



where, 



1 3 +3 / r 7(*)' 



with the following condition: 



3f) 
2 7 (t) 



(22) 



'5, 



m 



>izm -i e 



- E 



2^ 1 v^(t)5,(t); 2^5,(t) 

The semiclassical equation (17) can be rewritten as follows: 
5i(*)5 2 (t)5 3 (t) 



(23) 



8tt 



Z ' v Si(t) 5 2 (t) ^ 5 2 (t) S 3 (t) ^ Si(t) S 3 (t)>/ L 



f sinh 2 r + - 



E 



'Si(t)s s (ty 



4 4 fe v^(*)5,-(*). 



+iE 



WA 7,1/7 



4 ^ VS(*)/7 4V7 



- + - - +m' + 7 z 



+ (cosh r sinh r cos(c/? + 27^)) 



x 



E 



E 



7,1/7 



— I 7 I I + m* 



4 ^A^) Sj(t) J 4^V^(*)/7 4\7 
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(24) 



The above equation can be solved peturbatively. Starting from the approximation ansatz 
<Sio(t) = Slot™ 1 , S 20 (t) = S 20 t n2 , S 30 (t) = S 30 t n3 , and 7 o(i) = m , we obtain the next order 
perturbative solution for «Si, 

Stt™ I"/ _ 1\ / V 3 

S U (t) 



<S 20 S 30 (nin 2 + n 2 n 3 + nin 3 ) 
+ cosh r sinh r cos(c/? + 2mt) 



i=j=l n i n j 

8mH 2 



2— U2— nz 



(25) 
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Similarly, the next order perturbation solution for S 2 and <S 3 are obtained as 



S 22 (t) = 



8nm 



SioS 3 o [n 1 n 2 + n 2 n 3 + nin 3 ) 
+ cosh r sinh r cos(tp + 2mt) 



sinh 2 r H — 
2 



8m 2 t 2 



- 1 



and 



8m 2 t 2 
£2-m-n3 



(26) 



5 33 (t) = 



87rm 



5i 5 20 (nin 2 + n 2 n 3 + nin 3 ) 
+ cosh r sinh r cos((^ + 2mt) 



( sinh 2 



1 

r H — 
2 



1 + 



Si 



8m 2 t 2 



- 1 



8m 2 t 2 

2— Mi- ri2 



(27) 



Where «Sn means the next order perturbation solution for the scale factor in the x direction 
and the same hold for <S 2 i and <S 31 ; respectively in the y and z directions. 



From (25), (26) and (27) it follow that 



'11 



t 



2— r*2— 



'22 



t 



2— m— n3 



, S 33 ~ t 



2 —n\—ri2 



(28) 



Next goal is to examine the particle creation the anisotropically evolving universe described 
through the Bianchi type I metric. For this, consider the Fock space which has a one 
parameter dependence on the cosmological time t. Then the number of particles at a later 
time t created from the vacuum at the initial time t is given by 



N (t,t ) = (0,<f>,t \N(t)\0,4>,t ), 



(29) 



where N(t) = at a. Thus using (11), the vacuum expectation value of the right hand side of 
(29) becomes 

(N(t)) = (S 1 S 2 S 3 ) 2 4>4>*(4> 2 ) + #*(tt 2 ) - S&SzWiW) - S 1 S 2 S 3 ^*(^). (30) 
Therefore 

N (t,t ) = (^5 2l S 3 ) 2 |0(t)0(to)-0W0(to)| 2 - (31) 

The number of particle created in the vacuum states can be obtained by using the petur- 
bative solution in the limit mt Q , mt > 1, in (30), therefore, 



Wo(Mo) = 



SiS 2 S 3 



4 7(*)t(*o) S10S20S30 



SiSjn 



A CO C_ C. O 2.^ 



4 i=j= i \SiSj SioSjoJ 2 . =1 



J'O 
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1 3 <? 

+^E 



i(t) ] i(to)) 



+ 



1 (lit) 7(t )V 2 

2 7(W 



, If go fig) 7fa)' 
2S5 !0 U) 7(*o). 



+ 7W 2 -7(^o)' 



ni + n 2 + n 3 \ 2 /t-t \ 2 /f\ ni+n2+ " 3 
4^ J V tt ) \T , 



(32) 



By similar procedure, one can compute the particle creation for the quantized inflaton 
in squeezed vacuum states also. For this, consider the quantized inflaton in the squeezed 
vacuum state formalism. Then the expectation values of 7r 2 , <fi 2 , ircf) and (f>7i can be computed 
in squeezed vacuum state by using (11), (14) and (16), and are respectively obtained as 



) sqv 


= (SiS 2 S 3 ) 2 [(2sinh 2 r + l) <p*(t )<p(t ) + sinhr 


coshr ( e - l >* 2 (t ) +e% 




(0 2 ) sqv 


= (2 sinh 2 


r + l) (f>*(t )(f){t ) + sinhr cosh r (eT^ 


0* 2 (t o )+e^ 2 (t o )) 




{^4>)sqv 


= S1S2S3 


sinh 2 r<f)*(t Q )(j)(t ) + cosh 2 r<j)(t Q )(j)*(t ) + 






sinh r cosh r (e~ i,f <j>* (t ) (f>(t ) + e^0(t o )0(*o))" 




(33) 


{4>n) sqv 


— S1S2S3 


sinh 2 r0*(i o )</>(*o) + cosh 2 r0(£ o )0*(t O ) + 






sinh r cosh r (V^(t )<f (to) + e i{p <j>{t )j>{t Q ))~ 







Applying the above results (33) in (30), the number of particle created in squeezed vacuum 
in the limit mt , mt > 1, is obtained as 



N sqv (t, to) 



1 



S1S2S3 



47(^)7(^0) S l0 S 2 oS 30 

1 w) 

4 [ti Si(t) + 7 (t) 



( 2rfr+i )(Ksti + tl) 2+ ^ )+f(i) 

+ sinhr coshr (V^" 2 i*^)^) - e Hv-* S i(to)*o) 



-sinh 2 rfi (fM + M 

i v 2^^(t) + 7 w; vfeM*) 1 7(*b); 



E 



27(*)7(*o) 



(34) 



cosh 2 r(i (y*® + M 



El€ + ^)+2-,(*h(/o) 



^^o(t) 7(^0), 



sinh r cosh r (e"* + e^" 2 /**>)*<>)) ( - E lr 



2^(i) 7 (t), 



f 7(*o)Y 
vSSo(t) 7(to)| 
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Again using the peturbative solution , we get 



N sqv (t,t ) ~ 



{n 1 + n 2 + n 3 ) 2 /t\" 1+ " 2+ " 3 



4m 2 



+4m 2 sinh 2 r 



to 



(2 sinh 2 r + l) 



2t t 



sinh r cosh r cos (/? 



(35) 



Therefor the total number of particle created for the whole range of squeezing parameter 
and squeezing angle can be written as 



("Y'max f^pmax 

N sqv -tot = N (t, t ) + / N sqv _ r dr + / N sqv _ v d(p, 

Jo Jo 



(36) 



where 



N. 



sqv—r 



4 sinh r cosh r 



(ni + n 2 + n 3 y 
4m 2 

(sinh 2 r + cosh 2 rj cos </? 
2^ 



t-to 



2t t 

■j- \ rti+«,2+n3 

T 



2m 2 



(37) 



and 



AT. 



sqv — ip 



(ni + n 2 + n 3 y 
Am 2 



sinh r cosh r cos ip 
2M n 



-j- \ ni+n2+n:j 

T , 



(38) 



4. Conclusions 

In this paper we have examined the behavior a homogeneous and massive scalar (inflaton) 
field minimally coupled to the gravity in Bianchi type -I model of the universe, in the frame 
work of semicalssical theory of gravity. The inflaton is represented in squeezed vacuum state 
formalism of quantum optics and hence the approximate leading solution to the semiclassical 
Einstein equation is found. The next order solution for each scale factor in their respective 
direction show that each scale factor in each direction dependent on power law of expansion. 
Further more the solutions show that evolution of scale factors are mutually correlated. 
When 77.1 = n 2 = n 3 = n, then the corresponding solution reduces to isotropic model and is 
consistent with the result obtained as in the ref. 23. Form anisotropic to isotropic transition 
a damping mechanism is required. One of the efficient damping mechanisms could be due to 
the particle creation in anisotropic models. We have also examined the nonclassical particle 
creation in Bianchi type -I cosmological model by representing the inflaton in squeezed 
vacuum state formalism. The present study can account for the nonclassical particle creation 
and power law of expansion of the scale factors in Bianchi type -I universe, for a homogeneous 
and massive scalar field minimally coupled to the gravityin the frame work of semiclassical 
theory of gravity. 
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